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SECTION I 
INTRODUCTION 


This report describes a computer program to generate fatigue spectrum 
loading sequences. The report contains all information needed for the 
practical use of the program. 


The program is specifically tailored for the development of random 
cycle-by-cycle, flight-by-flight loading sequences typical of aircraft 
Structures. However, it can be used to generate a spectrum of any type for 
any structure. 


The random sequence of cycles and flights is produced by a random number 
generator. The basic loads occurrences spectrum is calculated through the 
use of the computer program A6PA (program 16PA in Reference 1). The descrip- 
tion of this program is not repeated here, but is to be found in Reference 1, 
with the exception of changes described in Appendix A. 


The program is written in Fortran IV language for CDC CYBER 74 computer 
or its equivalent. Douglas Aircraft Company program identification number 
is A6PD. The program was extensively used to generate loads spectra sequences 
in the Air Force study "Effect of Transport/Bomber Spectrum on Crack 
Propagation," Reference 2. 


SECTION II 
PROBLEM STATEMENT 


Structures subjected to periodic or non-periodic fluctuating loadings 
experience gradual degradation in its load carrying capability through forma- 
tion and propagation of fatigue cracks. Often, this type of loading is 
referred to as "fatigue loading". Fatigue loadings can be divided into two 
basic categories: 


1. Constant amplitude loading - fatigue loading in which all of the peak 
loads (maximum loads) are equal and all of the valley loads (minimum 
loads) are equal, 


2. Spectrum loading - fatigue loading in which all of the peak loads are 
not equal and/or all of the valley loads are not equal; this loading 
is also known as variable amplitude or irregular loading. 
Fatigue loading basic terms and examples of spectrum loadings are illustrated 
in Figures 1 and 2. 


Most fatigue loadings experienced by aircraft structures are of the 
Spectrum type. The spectrum features are a function of aircraft type and 
utilization, loading environments and type of structure. The sequence of 
loading is usually of random nature and is a function of loading environments 
and sequence of missions and flights. 


The basic loads spectra contents are usually supplied in the form of 
loads exceedances or occurrences information or other statistical data which 
do not supply any information about the sequence of these loads. An exceed- 
ances/occurrences spectrum specifies the number of times that the various 
loads in the spectrum are exceeded or occur in a specified period of time. 

An example of such information for Air Force aircraft is given in Reference 3. 


In view of the fact that the fatigue failure process in crack initiation 
and propagation stages is loading sequence dependent, means are needed to 
generate representative fatigue loading sequences which can be used in design, 
analysis and verification testing. The computer program described herein 
Provides the means of generating a random sequence spectrum from given exceed- 
ances/occurrences spectra, including the generation of flight random sequence 
from a specified mission-flight mix. 
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Figure 1. Fatigue Loading Basic Terms 


Figure 2. Examples of Spectrum Loadings. 


SECTION ITI 
PROGRAM OUTLINE 


The program name is Spectrum Loading Sequence Generation (or RANDOM) and 
its identification is A6PD. It consists of program A6PA, described in Refer- 
ence 1, and Appendix A, which generates the occurrences spectrum, and the 
additions to generate and edit random cycle-by-cycle loading sequences. A 
general outline of the program is given by Figure 3. A detail step-by-step 
explanation of the program, beyond the A6PA part, is given below. 


1. The basic information extracted from A6PA, to generate the loading 
sequence, are the segment-by-segment occurrences spectra (MIN,MAX,n). 
However, program A6PD contains the entire A6PA program and executes all 
its features. 


2. The occurrences spectrum generated by A6PA must be defined by the 
following A6PD input: 


NFT = Number of flight types. A flight type can be a realistic flight 
or any group of cycles to be identified as a "flight". (Card 2) 

NF = Number of flights for each flight type. (Card 4) 

NS = Number of A6PA segments in each flight type. (Card 5) 


3. All A6PA segments in each flight type must be identified by input flags 
Fl and F2. These flags are integers and must start with 1 and 
increase consecutively. Enter F1l=F2=0 for any A6PA segment to be 
excluded. If the cycles (n) of several A6PA segments with the same 
MIN,MAX values are to be summed, enter the same Fl flag for these seg- 
ments. If several A6PA segments with different MIN,MAX values are to be 
combined, enter the same F2 flag for those segments. A6PA segments are 
summed or combined, because in the loading sequence generation it is 
sometimes more realistic to treat several A6PA segments as one. Such 
situations might exist where several A6PA segments represent loadings in 
the same flight interval, such as maneuvers and gusts, but it is desired 
to have the cycles distributed within one F2 segment. In the sequence 
generation, the segment sequence is the F2 sequence. Example of com- 
bining A6PA segments under Fl and F2 flags is given on the following 
page. 


15,000 


10,000 15,000 


A6PA segment occurrences spectra are formed into a 
(MIN ,MAX ,n ,n/FLT) 


array according to flight type and F2 segment. All cycles n are 
rounded off to whole numbers and n<.5 are set to zero. 


If x(n/FIt) for a given F2 segment is not an integer, a ‘fictitious’ 
load level 
(MINe = MAX¢, ng) 


is added to the segment array to make £(n/FIt) for that F2 segment 

an integer, where, 

MINe = MAXe = (MIN] + MAX])/2 if © (MIN] +MAX])/2 = (MIN2 + MAXg)/2 
= MIN} if MIN] = MAX 
= MAX) if MAX, = MAXo 
= (MIN]+MAX) )/2 if none above. 


Next, selection must be made as to the type of flight sequence desired. 
Input IFRS and IAFS (card 2) and NFRS (card 12) if IFRS #0. The 
flight sequence may be: 


a) Random. IFRS=IAFS=0. The sequence is generated using a pseudo random 
number generator, see subroutine RANIC, page 88. The draw from 
the flight array is made without replacement. The standard seed 


number for the flight sequence random number generation is 11111 and 
is obtained by entering KF=0 in card 2. For a different random sequence 
enter 2051 < KF < 4194304 . 


b) Alternate Sequence. Sequence according to the largest peak per 
flight. IFRS = 0 
IAFS = 1 , LO-HI sequence 
= 2, HI-LO sequence 
= 3 , LO-HI-LO sequence . 


c) Specific sequence. IFRS = N , number of entries of ‘flight type 
and number of flights' sets entered in card 12 (NFRS). Also, 
enter IAFS=0. 


The spectrum loading sequence is generated by defining a sequence of 
valleys and peaks. Valey is MIN and peak is MAX as taken from the 
A6PA (MIN,MAX,n) array. A (valley, peak) pair is counted as a cycle, 
with valley preceding the peak. Groups of cycles are identified as a 
'flight' according to NS (number of A6PA segments in each flight type) 
and Fl and F2 (sequence of these segments). 


Choose one of two available valley/peak coupling methods: 


a) IRS=1. The valley/peak pairs are defined directly from the program 
A6PA output (MIN, MAX. , n.) array pairs, where 


Valley = MIN. and Peak = MAX. 
The sequence of cycles (valley/peak pairs) so defined is random 
and is generated using a pseudo random number generator, see 
Subroutine RANIC, page 


b) IRS=2. Valley = MIN and Peak = MAX are chosen individually and 
alternatingly in a random sequence from the program A6PA output 
(MIN,MAX,n) array. The random sequence is generated by a pseudo 
random number generator, see subroutine RANIC, page 88. 


In both of the above cycle sequencing coupling methods in using the random 
number generator, the draw from the (MIN,MAX,n) array is made without 
replacement. The standard seed number for this application is 12345 and 
is obtained by entering KC=0 in card 2. For a different random sequence 
enter 2051 < KC < 4194304 . 


The cycle (valley, peak) sequence established in the previous steps is 
edited in the following order by EDIT] and EDIT2 : 


EDIT 
All valleys and peaks which are not loading reversal points are 
eliminated. This effectively eliminates fictitious load levels 


(if they do not produce a loading reversal) as well as intermediate 
data points. 


EDIT2 
a) Range Truncation. Elimination of cycles as a function of cycle 
range (peak-valley) and R (valley/peak). If no range truncation 
is desired, NXY=0, card 2. If range truncation is to be per- 
formed, define the range truncation level over the desired R 
interval by XY input in card 9 


data points all cycles with (RANGE,R) 
values below this line 
eliminated 


—— ae 
ALOT 
YMMV 
and specify the number of data points (RANGE,R sets) entered as 
NXY in card 2. The program interpolates linearly between the 
data points and will take the closest RANGE entry if cycle R 


is outside the R interval entered. In the latter case an 
ERROR statement will be printed. 


RANGE 


There are two exceptions to this range truncation procedure: 


1. When a peak or valley, from the cycle to be eliminated, is 
higher or lower than the preceding and next peak or valley. 
In such cases the subject peak or valley is retained, but 
their companion valley or peak, as well as the preceding peak 
or next valley, are eliminated. 


| retain 
final cycle po retain — 
eliminate—mp yi; ____.__ candidate cycle aaa mae grees 
‘| for elimination vA 2litaxes 


10. 


Lee 


2. When a cycle to be eliminated (valley., peak, ) is followed by 
a cycle with valley... > valley, and peak... > peak, 
In such case peak, and valley.) are eliminated and the 
resulting cycle becomes valley, peak, ,, ‘ 


——-retain 


candidate cycle p final cycle 


for elimination Y eliminate 


b) Cycle elimination if the peak is larger than a specified value 
ELIMP, card 3. 


c) Clipping, according to specified peak (CLP) or valley (CLIV) 
magnitude, card 3. Given a cycle (valley, , peak.), 


1T peak. > CLP , set peak, = CLP 

if peak, > CLP and valley, > CLP , eliminate cycle 
1T valley, < CLIV , set valley, = CLIV 

1T valley, < CLIV and peak, < CLIV, eliminate cycle. 


d) Multiplication of the complete valley, peak sequence spectrum by 
a constant, FACTOR in card 3. Enter FACTOR=1 when multiplica- 
tion is not desired. Multiplication will not be performed if an 
alternate (IAFS=1, 2 or 3) flight sequence is used. 


If, through preceding editing, all cycles in a flight are eliminated, the 
flight will be counted in the total, but it will not be included in the 
output printout nor magnetic tape. An ERROR statement will be printed. 


The basic output of program A6PD is a flight-by-flight loading sequence 
of valleys and peaks. For detail description of the output see the 
section on OUTPUT, page 101. 


The program will terminate according to input IPTF (card 2): either when 
all flights (specified by NF, card 4) are processed or after N_ specified 
number of flights. 
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Variations in a spectrum can be accomplished through the various features 
provided in the A6PD input data as well as by changing the basic occurrences 
spectrum from A6PA through A6PA input data. Some of the variations that can 
be accomplished through the A6PD input data are: 


1. Combining and/or elimination of A6PA segments. 

2. Different Valley/Peak coupling. 

3. Different flight sequences. 

4 Elimination of cycles as a function of (RANGE, R), valley or peak values. 
5. Clipping of loads below or above a specified value. 

6. Multiplication of all loads by a constant. 


7. Variation in spectrum length in terms of number of flights. 
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[AGPA | 
| INPUT AGPA 


LDATA _J 
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SPECIFIED BY A6PD DATA 


SEQUENCE EDITING: 
@ ELIMINATION OF INTERMEDIATE DATA POINTS BETWEEN REVERSALS 
@ TRUNCATION (ELIMINATION OF CYCLES AS A FUNCTION OF 
(RANGE, R), PEAK OR VALLEY) 

e CLIPPING (LOADS BELOW OR ABOVE A SPECIFIED VALUE SET 
EQUAL TO THAT VALUE) 

@ MULTIPLICATION OF ALL LOADS BY A CONSTANT 


OUTPUT 
PRINTOUT FOR A 
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Figure 3. Spectrum Loading Sequence Generation Program (A6PD) General Outline. 
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SECTION IV 


DEFINITION OF SYMBOLS 


The A6PD input data symbols and definitions are given in Section V. 


Following are the additional symbols used in the computer program. 


A6PA symbols and definitions see Reference 1 and Appendix A. 


A (or N) 


ACY 


AFS 


HPEAK 
ICASE 


IDPEAK 


IERR 


IFLAG 


IF 
IF2 (or IS2) 


JJ 


IMM 


The dynamic working array. Throughout the program core 
allocation is done through implied equivalence to this array. 


An array of cycles read in fron A6PA. 


The reordered array of highest peak per flight (alternate 
flight sequence). 


An array of maximum stresses read in fron A6PA. 

An array of minimum stresses read in fron A6PA. 

An array of cycles used during combining of segments. 
An array of the highest peak per flight for printout. 
A6PA case number. 


An array of flight number identification for the highest 
peak per flight printout. 


A running sum of the number of error messages written. 


An array used in reordering the flight sequence to indicate 
when a flight has already been chosen. 


An array of Fl segments. 
An array of F2 segments. 


An array used for debug that sums the number of cycles 
dropped out for the various editing cases. 


A cycle counter used during edit type 1 and 2. 
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For program 


IRR A6PA reference run number, or, an array used for the 
cumulative totals in the spectrum summation table. 


ISEG A6PA segment number. 


IS] The beginning segment number for printout of flight type, 
reference run, case and segments, or, the beginning segment 
subscript of the IFl and IF2 arrays. 


ISz The ending segment number for printout of flight type, 
reference, case and segments, or, the ending segment 
subscript of the IF] and IF2 arrays. 


JTN An array of the flight types corresponding to a flight number. 
KCY An array of cycles after combining of segments, used 

when IRS = 2. 
LEFT The amount left of NSIZE usually calculated before a 


subroutine call after array sizes have been set. 


LREC An array of reordered flight numbers for the alternate 
flight sequence corresponding to the AFS array. 


MCY An array of the number of cycles per flight. 
NCY An array of cycles after combining of segments. 
NFF An array containing the sum of the number of times a 


flight type was picked and the sum of cycles for each 
flight type. 


NPG A page count. It is written at the top of each new 
output page. 


NSIZE The dimension of the "A" array. 
NST The total number of segments. 
NTF The total number of flights. 
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PMAX 


PR 


RANGE 


SMAX 


SMIN 


SMM 


An array of the highest peak per flight used for the 
alternate flight sequence. 


An array used for the spectrum summation print of peak 
WSs K« 


Cycle loading range = maximum stress (peak) - minimum 
stress (valley). 


Stress ratio = minimum stress (valley) / maximum stress (peak). 


An array used for the spectrum summation print of 
range vs. R. 


An array of maximum stresses used during and after combining 
of segments. 


An array of minimum stresses used during and after combining 
of segments. 


The final array of valley (minimum) and peak (maximum) 
stresses sequence. 
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SECTION V 
INPUT DATA 


The input data is divided into two groups: 


1. Program A6PA data required to calculate the exceedances/occurrences 
spectra. This input data and data sheets are described in Reference 1 
and Appendix A. 


2. Sequence generation data (henceforth called "A6PD data") to generate the 
sequence of cycles and flights and to perform other manipulations in the 
A6PD program. The description of this input data and the sequence in 
which it is input follows. 


CARD 1 - TITLE = Enter job description. Use only one card. 
CARD 2 - NFT IRS IFRS IUIL NPI NRAN KVP IFI NXY IPFS NPSS IPTF IAFS MAXHP KF KC. 


NFT = Number of flight types. 
IRS = Valley/Peak coupling format. 
= 1 3; program A6PA pairs 
= 2 3; individual (random) 
IFRS = Flight sequence. 
= 0 3; random. Also when IAFS # 0. 
= N 3; number of ‘flight type and number of flights' sets 
when special sequence is specified. Enter the ‘flight 
type and number of flights' data on card(s) following 
Fl and F2 data (card 12). 


IUIL = 1; Fortran unit number 
ea = Number™ of peak levels in spectrum summation 
TRAN a Wten” a _—e 
NRAN = Number of range levels in spectrum summation 
IVP = Number” of R levels i 
KVP = umber 0 evels in spectrum summation 


* Enter the number as a negative number if the peak, range 
or R_ values are input directly in cards 6, 7 or 8. 
A positive number means that only the extreme values are 
entered in cards 6, 7 or 8 and the levels are calculated 
at equal intervals using the NPI, NRAN or KVP inputs. 
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IFI 


NXY 


IPFS 


NPSS 


[FTP 


IAFS 


MAXHP 


IKF 
KF 


IKC 
KC 


Fortran unit number. 

0 ; spectrum output tape not to be written. 

3 ; spectrum output tape to be written. 

Number of data points (R, range sets) in R vs range input 
data (card 9) if range truncation is to be performed. 
NXY=0 if no range truncation is to be performed (no 
data is to be input in card 9). 

Print loading sequence? 

-1 3 no. 

0; yes, all flights. 

N ; yes, N flights. 

Number of spectrum summations to be printed. In card 10 
indicate the flight number(s) for spectrum summation(s). 
Flight number at which the program is to be terminated. 
0 ; normal, equals to total number of flights shown on 
card 4. 

N ; stop after N flights (less than the total shown on 
card 4). 

Alternate flight sequence. 

0; random. Also, when IFRS=0 or N. 

1 ; LO-HI (on the basis of the largest peak load 

per flight). 

23; HI-LO 

3; LO-HI-LO 

Number of largest peaks per flight to be printed 


Seed number in the random number generator for 

sequencing flights. 

0 ; will default to the number 11111 which is considered 
to be the standard number for this program. 


Seed number in the random number generator for 

sequencing cycles. 

0 ; will default to the number 12345 which is considered 
to be the standard number for this program. 
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CARD 3 - 


CARD 4 - 


CARD 5 - 


CARD 6 - 


CARD 7 - 


CARD 8 - 


CARD 9 - 


CARD 10 - ISS 


CARD 11 - F1,F2 


CLIP CLIV FACTOR ELIMP 


CLIP = Peak clipping value. 

CLIV = Valley clipping value 

FACTOR = Multiplication factor 

ELIMP = Peak value for elimination of cycles with peaks 


above this value. 
(Input very large numbers for CLIP, CLIV or ELIMP to prevent 
clipping or cycle elimination. ) 


NF = Number of flights for each flight type that the A6PA 
data represents. Enter in increasing flight type order. 


NS = Number of A6PA segments in each flight type. 
Enter in increasing flight type order. 


PI = Peak values for spectrum summation. 
Enter all values if NPI in Card 1 is entered negative. 
Enter only the two extreme values if NPI in Card |] 
is entered positive. 


RAN = Range values for spectrum summation. 
Enter all values if NRAN in Card 1 is entered negative. 
Enter only the two extreme values if NRAN in Card 1 
is entered positive. 


VP = R values for spectrum summations. 
Enter all values if KVP in Card 1 is entered negative. 
Enter only the two extreme values if KVP in Card 1 is 
entered positive. 


XY = Rvs range data (R, range sets equal to the number of 
sets specified by NXY in Card 1). Enter data in R 
increasing order. 


Flight number for spectrum summation. If making more 
than one entry, enter in increasing order. 


Flags defining Fl and F2 segments. Enter Fl and then F2 
flags for each flight type. Enter data in flight type 
increasing order. Start each set of Fl and F2 data for 
each flight on a new card. 


2) 


CARD 12 - NFRS = Flight sequence when IFRS # 0 in Card 2. 


Enter flight type and number of flights sets in the 
desired sequence where the number of sets is equal 
to IFRS=N. 


General Notes about A6PD Input Data: 


\ 


The data consists of a string of values separated by one or more blanks 
or a comma. Blanks are not allowed within a constant. A decimal point 
omitted from a real constant is assumed to occur to the right of the 
rightmost digit of the mantissa. 


Except for CARD 1 data, more than one card can be used to enter a 
particular data group identified as CARD X data in the preceding 
description. If more than one card is used for a given type of data, 
consider the cards as 6.1, 6.2, 6.3 or 11.1, 11.2, 11.3, etc. Input 
cards cannot have line sequence numbers. 


Enter the data in the sequence indicated in the preceding writeup. 
This sequence is summarized on the following page in the form of an 
'A6PD INPUT DATA FORMAT’. 


An example of A6PD output input data printout for a sample case is 


shown on pages 208 and 209 . 
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SECTION VI 
PROGRAM LISTING 


The following pages contain the program source language listing generated 
by the Fortran Extended (FTN) compiler. 


Program A6PA (Ref. 1 and Appendix A) listing, as adapted to the CDC 
computer and A6PD program, is also presented here. 


Figure 4 shows the interaction of all the subroutines. The listings 
of these subroutines, including additional descriptions, are to be found on 
the following pages: 


Ne DESCRIPTION PAGE 


eat Calls subroutines MMAIN (program A6PA) and RDMAIN 28 
Program A6PA subroutines: 
MMAIN A6PA main program (includes BLOCK DATA) 29 
NPUT1A Reads A6PA input data. 4] 
SPECSM Performs segment spectra summation. 48 
PRINT Prints A6PA input data. 90 
OMEVAR One variable interpolation. 97 
TWwOVIN Two variable interpolation. 59 
Spectrum Loading Seqiience Generatton Subroutines: 
RDMAIN Reads and prints title and other A6PD input data. 63 
INPUTF Reads A6PD input data. 65 
INF1F2 Reads and prints remainder of A6PD input data. 66 
NEWPG Prints output page heading and page number. 70 
INMMN Sets up core storage required for AMMN. 7| 
AMMN Combining of segment spectra using flags Fl and F2. 72 
OPENMS Opens mass storage file. 76 
GENFL Generation of flight sequence: random or specified. a 
GENAFS Generation of alternate flight sequences on the 80 
basis of highest peak per flight. 
GENCY Generation and editing of the random cycle sequence, 83 
DISTRD Selection of flight or cycle number through the 87 
use of the random number generator, RANIC. 
RANIC Random number generator. 88 
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Ne DESCRIPTION PAGE 


REED Reads the sequence of valleys and peaks for one 89 
flight. 

READMS Transmits data from mass storage to central memory. 90 
REDIT1 Repeats EDIT1 type editing found in GENCY. 9] 
SPSUM Performs spectrum summations. 92 
PRNTSS Prints spectrum summations. 93 
ERROR Prints error messages. (See page 103 for a listing oF 
of the messages. ) 

RITE Writes spectrum for one flight into temporary 98 
storage. 

WRITMS Transmits data from central memory to mass storage. 99 
WTAPE Writes spectrum on output tape. 100 
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SECTION VII 
OUTPUT 


Program A6PD output (excluding A6PA part which is described in Reference ] 
and Appendix A) consists of the following printout: 


1. Input data. This also includes the relationships between flight type 
number and program A6PA RR, case and segments. 


2. (MIN,MAX,n,n/FLT) arrays for all flight types and F2 segments. 


3. Flight-by-flight valley, peak loads sequence for a specified number of 
flights (IPFS input, card 2). 


4.  Sepectrum summations 
(RANGE, R) vs (n,=n) 
and (PEAK, R) vs (n,n) 
The input required to obtain this output is: NPI, NRAN, KVP, NPSS 
(card 2), PI (card 6), RAN (card 7), VP (card 8) and ISS (card 10). 
A number of summations, after different number of flights, can be 
obtained by appropriate NPSS and ISS inputs. 


9. Number of flights and cycles for each flight type in the spectrum 
generated. 


6. The magnitude of the N number of largest peaks per flight in the 
spectrum generated, together with the flight number in which they occur. 
Input N = MAXHP, card 2. 


The complete valley, peak loads sequence may be written out on magnetic 
tape for further use in fatigue or crack propagation analysis or testing. 
Input IFI=3 (card 2) if spectrum on tape is desired. The tape contains, first, 
the total number of flights, and then, flight number, and the number of cycles 
and the valley, peak loads sequence in that flight, for all the flights. 


For sample output printout see pages 208 thru 235 . 
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SECTION VITI 
ERROR CONDITIONS 


There are seven special error conditions written for program A6PD. Each 
error is written out in the following format: 

ERROR N IS TYPE M AND INVOLVES J, K, L, 

where, N = error message number; the program will stop execution 


when N exceeds 900 and will print the completion 


code message STOP 900. 


M = error type number. 


J Keb special description for each error type. 


The error conditions and probable steps to be taken to correct the 
condition, are: 


TYPE 
1 


DESCRIPTION 


More core is required for this job. The dimension of the A array 


in Subroutine RDMAIN and the value of NSIZE which is the 
dimension of A must be increased. 


J = amount of additional core needed 
K = the NSIZE for this run 
NOTE: Add 2*MAX# of cycles to "J" also. 


The numbering of the F2 segments is incorrect. 


overall segment number 
number of F2 segments up to this point 
the number of Fl segments up to this point. 


L 


More core is required for this job. (Message written in 


Subroutine INMMN.) The dimension of the A array in Subroutine 


RDMAIN and the value of NSIZE, which is the same as the 
dimension of A, must be increased. 
J = the amount of core needed at this point 


K NSIZE specified for this run 
M flight type number. 
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More core is required for this job. The program is terminated, 

and a completion code message STOP 9001 jis written. (The messages 
are written in Subroutine GENCY.) The dimension of the A array 
in Subroutine RDMAIN and the value of NSIZE, which is the same 

as the dimension of A, must be increased. Look at the indicated 
flight type. 


J = number of cycles allowed for in this run 
K = amount of total core specified for this (NSIZE) 
L = flight type that exceeds the core requirements 


The Valley to Peak Ratio of a cycle exceeds the maximum Valley to 
Peak Ratio input. The cycle is not eliminated. If this is not 
acceptable, correct the input for the Valley to Peak Ratio 

versus Range elimination curve. (This error will not be printed 
beyond five times in one run.) 


J = flight number 
K = cycle number 
L = 1] first end of curve 


2 second end of curve 


A flight has been edited until it has zero cycles left. The 
flight is not written on tape. 


J = flight number 
K = flight type 
L = undefined. 


In Fl segment (MIN,MAX,n) array cycles (n) for all load levels 
are zero. Either enter cycles for those A6PA segments making up 
the Fl segment, or use Fl and F2 flags to eliminate them. 


J = Fl segment number 
K = total number of segments in the flight type 
L = total number of A6PA peaks, valleys and cycles for this segment. 
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SECTION IX 
KEYPUNCH AND DECK SETUP INSTRUCTIONS 


Keypunch Instructions 


The following instructions apply to A6PD input data only. See Reference 
1 for A6PA input data instructions. 


1. Use standard keypunch procedures and a standard 80 column data card. 
2. All 80 columns of the data card can be used to enter the data. 


3. Except in card 1, which is the job description (title) card, see 
Section y, the data entries must be separated by one or more blanks or a comma. 


4. Data cards must be stacked in order by the user, as defined in 
Section V. 


Deck Setup 


The program and input data deck setup, where the input data is entered 
from a punched card deck, is shown in Figure 5. The input data represents 
information to generate one spectrum loading sequence. An alternate operation 
procedure, where the input data is read in through CDC utility UPDATE, is 
given in Appendix B. Job control card content and arrangement varies depending 
on the particular computer facility and procedures used by the user. Douglas 
Aircraft Company computer facility and procedures require the following job 
control cards to run on the CDC computer running under KRONOS: 


SER, ALQ1. 
A6PD, MFL = 155000, Ty | 4 
ACCOUNT, 2 


CID, Name, etc. 

FTN, I=INPUT ,L=9. 

C@PYBR( INPUT ,DATA1), 

REWIND ,DATA1. 

REQUEST , TAPE3 ,NT ,LB=KU ,D=1600, VSN=RESERVE ,T 
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LGQ,DATAI. 


Where, 


In Card 2 


In Card 7 


MFL 
- 


1/0 


T 


core size required for the program = 155000. 


estimate of CP time (seconds) to run a job. 
(Allow approx. 50 seconds per 10° cycles of 
the A6PA spectrum. ) 


estimate of input/output processing time (seconds). 
(This time varies significantly according to the 
scope of the job. Allow at least 15 seconds per 
10° cycles of the A6PA spectrum. ) 


estimate of line printout (in thousands). 
(Printout will depend on the size of the job, the 
extent of the optional A6PA output and the length 
of final spectrum sequence printout. ) 


Number of days that the output tape is to be 
reserved (saved) in the tape library. 
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re 4 


A6PD INPUT 
, DATA 
AGPA INPUT 

DATA 

ND OF RECORD 

——___x > 
AGPD PROGRAM 
END OF RECORD 
JOB CONTROL CARDS 


Figure 5. Program A6PD Deck Setup. 
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SECTION X 
TEST CASE 
This section presents a test case to check and illustrate the capabilities 
of program A6PD. The case represents spectrum BS1 from Reference 2. 


The spectrum represents two flights of a STOL transport basic and 
alternate employment mission: 


PROGRAM A6PA NUMBER OF FLT.HRS. FLIGHT A6PD 
FLIGHT | FLTS. * LDGS. HOURS FLIGHT TYPE 


1,099.8 
] cy 4 


TOTAL : 


Program A6PA segments and Fl and F2 flags for A6PD segment definition 
are described in Table 1. The spectrum represents wing root lower surface 
stresses. For more background information about these flights, see 
Reference 2. 


The A6PA and A6PD input data is presented on standard load sheets, pages 
112 through 165 . The S-N input data (data load sheet V-1.1, page 128) 
in program A6PA is fictitious so as not to produce any fatigue damage calcula- 
tion since the objective of this case was only to generate a spectrum loading 
sequence. 


The output is presented on the following pages: 


Pages 
A6PA, - Input data 112-164 
- Qutput 166-207 

A6PD, - Input data 165 
- Qutput 208-235 


The following output is not shown: 


1. A6PA individual segment and summary damage. The damage is zero because 
of the fictitious S-N data. 


2. A6PA, RR 15, Case 2 output. The data is very similar to Case 1, 
except for differences shown in Table 1 and higher stresses. 
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Note that, although the spectrum generated represents 2,528 flights, 
only 50 flights of the loading sequence printout were requested. 


Time to run this test case on a CDC computer was: 


I0 = 25.0 SC. 
CP = 100.3 sec. 
MRU = 1.8 
SRU = 12.5 
Lines printed = 4,181 
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TABLE 1 
TEST CASE A6PA SEGMENT DESCRIPTION AND F1 AND F2 FLAGS 


Fit. 1A Fit. 1B 


A6PA 
Pre-Fit. Taxi 
Climb, FD* Gust 
Climb, FD, Manv. 
Climb, FD, Manv. 
5 
6 
7 
8 
9 
10 
ll 
12 


Cruise, Gust 
Cruise Manv. 
Cruise 


Descent, 
Descent, Gust 
Descent, Gust 
Descent, Gust 
Descent, Gust 
Descent, Gust 
Descent, 
Descent 


Descent ,FB,Gust 
Descent,FD ,Manv. 
FD 


Manv. 


Descent,FD,Gust 
Descent ,FD,Marv. 
FD 


4 
4 
4 
4 
4 
4 
4 
§ 
$ 
6 
7 
8 
8 
8 
8 
8 
0 
‘ 8 
: 9 
1.5-0 
1.5-0 
1.5-0 


Descent ,FD,Gust Not Used 
Descent ,FD,Manv. Not Used 
Descent,FD,Manv. Not Used 


Landing Impact 
Post-Flight Taxi 


SL 
SL 


pet tee | 
mm 

” 
— 
“NO 


* FD = Flaps Down 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page__ of. 


FORTRAN DATA LOAD SHEET 1-3 Jasons 
STANDARD DATA INPUT rae ———— 
Pete) CKEYPUNCH: DO NOT 
lool) 6 pa} SEYPUNCH: PUNCH BLANK 
nT; [1.6] DATA FIELDS. 


L1} = FLAG WHICH SPECIFIES WHETHER THE MULTIPLICATION FACTOR F IN THE CYCLIC 
LOAD EQUATIONS WILL BE USED. 
= 1, YES. (COMPLETE SHEET III-1.2) 
= 2, NO. (F = 1.0) 
IW3 = FLAG WHICH SPECIFIES WHETHER THE CONSTANT P IN THE CYCLIC LOAD EQUATIONS 
WILL BE USED. 
= 1, YES. (COMPLETE SHEET III-1.3) 
= 2, NO. (P = 0 
IW] = FLAG WHICH SPECIFIES WHETHER GAG CYCLE INFORMATION IS TO BE PRINTED. 
ae ONLY IF 14 # O. 
= 2, NO. 
IW2 = FLAG WHICH SPECIFIES WHETHER THE SEGMENT SPECTRUM AND DAMAGE INFORMATION 
IS TO BE PRINTED. 
= I, YES. 
= 2, NO. 
IW4 = FLAG WHICH SPECIFIES WHETHER SPECTRUM SUMMATION FOR DIFFERENT ENVIRON- 
MENTS IS TO BE PRINTED. 
1, YES. (COMPLETE SHEET II-7) 


2, NO. 

FLAG WHICH SPECIFIES WHETHER THE INPUT DATA WILL BE PRINTED. 
1, YES. 

es NO. 


IW5 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page___of 


STANDARD DATA INPUT 1 Date 
69 70 7) 73 «77 60 


KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 


M5 = FLAG WHICH SELECTS THE SPECTRUM INCREMENTAL LOAD (ay) INPUT FORMAT. 
A VALUE (1 TO 7) MUST BE ENTERED FOR EVERY SEGMENT BEING ANALYZED, 
EXCEPT FOR THOSE SEGMENTS WITH M3 = 10, 11 OR 12. 


= |] TO 64 dy TABLES 1 70 6G. 
= 7; dy CALCULATED BY EQUATION ay = dy, + dy), CT mt Wels FE he Be Suet 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page____—ioff 


FORTRAN DATA LOAD SHEET II-3 Prepared by | 
STANDARD DATA INPUT 1 Date 

69 T0 Ti 73: PF 60 

1 Slooo KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 


RR CASE PROG 


Cae ee al al cel a I I a 


eee ee ee te ee ee ne ee we ee ne. me ae 6 ee ee ee we oO 6 0 6 6 6 8 se ete we oa gw be 6 we we Oh we 
OO OCONEE, LICL CI ELT RO I Pe Ee ee 
ee Ce eee he ee Be hee be eee ee he ee ee a Oe 
MMI MRM > RR ee eke the eee ee See ee hee ee he ee he eae he 


ss 
Yiwr 
nif 
me) 
Mm 
wn 


Seqm 11... 6 eqm 2] | “Ol XX 
or peer “pal 27), 
— Hl Py EE 

TY poet 
[sl oll... tt [oc esl all. 
pat ual WC 26f.aall. 
196 PY) SC 
Re ree oC Derren © ee eee Peer ie: 
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meer perererer 6: ae Bee peeeere 
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Tati Da ll 3711, ll 
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tl IE : » 
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ISTRES = FLAG WHICH SPECIFIES WHETHER STRESS TABLES WILL BE USED OR NOT. 
= 0, NO. 


= 1 TO 14, YES; A VALUE, 1 TO 14, CORRESPONDS TO THE STRESS TABLE 
NUMBER TO BE USED. (STRESS TABLES MAY BE USED CNLY WHEN 
mS = | 109, OR 13 TO1LS,) 
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FORTRAN DATA LOAD SHEET 11-4 Prepared by 
STANDARD DATA INPUT 1 Date 


V 


— 

|g 

— 

a4 418 ON WIT 

me Ser ere eer De 

gl wae |e a go ey el er eee 

ieee errr Fb: es 

re th fs ocen se call OO |.85) idiot ly yas enh OOM 

ified tt -—36livanat tas 
aed ALIX) ©[ 37, 342d 

Asad XX «6st ass EWE Ls 

A Peel eal. ss sana Doe 

LOO 68 | wna ee ae yy wee Ls 


NIFLAG = FLAG WHICH SPECIFIES THE LOAD CYCLE FORMAT. 
»= 1, F ((Y. + sy/2) + dy/2] + P 
= 2, F [Y. tay) + P 
= 3, F Lite - ay/2) + Ay/2] + P 


A VALUE MUST BE INPUT FOR EVERY SEGMENT BEING UTILIZED, EXCEPT FOR 
SEGMENTS WITH M3 = 10, 11 OR 12. 
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FORTRAN DATA LOAD SHEET II-5 
STANDARD DATA INPUT 1 


69 70 Ti 7S 77 60 


KEYPUNCH: 


Page of 
Prepared by 
Date 


DO NOT PUNCH BLANK DATA FIELDS. 


IA = 


FLAG WHICH SPECIFIES THE S-N DATA TO BE USED. 


A VALUE CORRESPONDING TO THE S-N DATA TABLE NUMBER AND DATA FORMAT AS 
DESCRIBED BELOW MUST BE ENTERED FOR EVERY SEGMENT BEING ANALYZED. 


7? TO i2 
3 TO 18 


19 TO 24 
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Y-ARGUMENT 


a EE EE 


cis Swax/SuLT Smin/Smax? = 
Swax/SuLT SMEAN/ SULT 
Satt/uLT SMEAN/ SULT 


SmMIn/°ULT 


FATIGUE DAMAGE CALCULATION PROGRAM 16PA 


FORTRAN DATA LOAD SHEET III-1.1 
STANDARD DATA INPUT 1 


69 70 71 73 77 80 
1,510.00} ..6 7.4] 


RR CASE 


PROG 


ve pete ee oO 0 Bb be 8 0 6 0 0.0 6 0.0 Om 6 6.8.0 0 0 0 0 0.0 0.8 0.0.0 8.0.0.0,8 
ote rere abe et eb ar6 01s 510 6.6.0 Me 6.0.0 MD 0.4.6.0. 0.0 8&4 6 05.0 0.0.8 00050589 
oretere ses 60 8 0.6 68 6 6 4 0.6 8 0 8 
IMI OT a Pte RM i mM NI A Oa te ee) PC 
OPO Rt, SD Sethe at te Re ee a od CR, a 
Pe Bea tO at ht eta eR Te 


SP it OD SO 
ntetetete tote ene ete lete ete ete ene o Ge 00.6 Be 6 6. 4.0.0.8.0 0.0.4 0 + 00.9, 0,4,6.6,0.0,0,° tetete oe 8 0 0.0.6 0.8. 0.9.8 
oto ea oe Gb oie ee 6 4 be 0 0.6 6 0.0 0.4.6 6.6.8 0 6.0 60 4 808 


. 
os 


— 
=] 
Oo 
bes 


m— Pa | 
Jjyi[n 
nlalo 
s[al[u 


Lv.6 

1,7,6 
mM slaty 
mH phi Cyt 


=e 


ca 
— 
=] 
=~) 
~m |= 


AM = CONSTANT LOAD, Yoo IN THE CYCLIC LOAD EQUATION. 


KEYPUNCH: 


Page. sof 
Prepared by 
Date 


DO NOT PUNCH BLANK DATA FIELDS. 


quan | Loc 
fe 


Mm 


— | re 
~~ 
~~ 
nf > 


_— 
oJ 
~ 
Oo 


WwW 


3 


W |W W |W |W ia] ia) 
“Sw [OV & |W IM |— |O ~ Ww | Ph 


& iw WwW 
an) le 


: as 
(MAY BE MIN., MEAN, 


OR MAX. LOAD DEPENDING ON THE LOAD CYCLE FORMAT.) 
A VALUE MUST BE ENTERED FOR EVERY SEGMENT BEING ANALYZED, EXCEPT WHEN 


M3 = 10 TO 12. 
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FORTRAN DATA LOAD SHEET III-3.1 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 Ti 73 «77 60 
is[ooolfi.e pA] KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 
rua [100 
Segm 21] , 2,7,6 
: ma) 
— i Serres: se ee 
— y) Serer asc ee 
_—: Jepmeeeuet ils | Se 
ae ee Serre 
ee eo eee i se) 
me Pe) eee a 
ae Pe re eee | 
ee Pe ee 
a ee rere ee 
ee ae) 
13 Hl LLLP | oe 
ee eee 
mae) 


eT PP ee 


ay, = FIRST TERM IN THE INCREMENTAL LOAD EQUATION ay = ay, + ay,, (J - 1.0) 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M5 = 7, 
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FORTRAN DATA LOAD SHEET III-3.2 Prepared by 
STANDARD DATA INPUT 1 Date 


ayaa | rae a oe =— 

SS ee Pee es 
a Per: 100} | 2 | eer T) s 
| et AOE ete eee VE) boy y) eres ts 
ey 7-1 OO 0 Oe TET) | PPO 
mae ee Per i eae ee 1 eT 
mae ee eee i a ee | eer a 
ee ere Tie ee ee ere Li 
cent cant LED) (Pel ee i gee 
er YY! TH) (ae a 
a ii oe re ee is 
ou 6 iit hcl prererritis 
12 FETT C32] y32y7] yg a | OO | 
13 Ce) te a ee ee ee 
14 Wit ee eer 
| 15 Nigee ol eee is 
ET: nite hl eres 
_: LET. LPP rrr ire 
| 18 ee —* reer irs 
_ 19 Le | ree 
Le lll 100 (rer 1 on 


by," SECOND TERM IN THE INCREMENTAL LOAD EQUATION Ay = Ay, + Sy), (j-1.0). 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M5 = 7. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA 
FORTRAN DATA LOAD SHEET ITI-4.1 
STANDARD DATA INPUT } 


69 70 7) 73> FP 80 


CASE PROG 


watatd 
ae Pee Ea 
SSO OO oe 


IZG 0000 


104] 
| 1%6 
loz 
ae 
ch 
na 
Ty 
10,2. 
lon 


s|-|-l-Isisisisis 
NET TTP fo fe fis fo 


In INPUT. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 7 TO 9, OR 13TO1S, 


KEYPUNCH: 


Page of 
Prepared by 
Date 


DO NOT PUNCH BLANK DATA FIELDS. 


+| VALUE 


Gua 


men eRe TD ER ae COR Bt et I 
Se ae ee Ce Be ke ee 

‘ene e*e ee aM ea tete’ ter eme B's ee rn 

anes * ate Hu x 


7 ae Y woes EEE — 


oe a 

are) eee 

| 2a], 35,9] 11,8,0,6 

as) 360 lee 3, 

2s 

36,3} 18)8)3) 1 ay 
564 Ry 5 a wey | 


Sis SLO SIS [Sd 
L21N <=; N LIN | & 


tn AT Ay = 0 FOR THE FIRST TERM OF THE GENERAL OR GUST EQUATION SPECTRUM 
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FORTRAN DATA LOAD SHEET III-4.2 Prepared by 
STANDARD DATA INPUT 1 Date 


69 70 TI v3. 77 60 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RR 


case PROG 


+| VALUE + 
Soo 


tes 12 2 
PT Mew ge x eg BO 
3,7 


40,2 
4,0,3 


38 
39 


No, = xn AT Ay FOR THE SECOND TERM OF THE GENERAL OR GUST EQUATION SPECTRUM 
In INPUT. ENTER No, = 0 IF SECOND TERM IS NOT USED. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 7 TO 9, OR 13 TOUS. 
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FORTRAN DATA LOAD SHEET III-4.3 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 Ti 73 77 80 


KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 


4 2,8 


~o 


4 3,0 
iy A342 


No; = fn AT Ay = O FOR THE THIRD TERM OF GENERAL EQUATION SPECTRUM fn INPUT. 
ENTER No, = 0 IF THIRD TERM IS NOT USED. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 7. 


>> > ih PP PP PP 
MN Mm NX JO [PO 
Oy im JP fw 


i 
| 
Ww 
t~ 


> 
Ww 


AN 
Ww 
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FORTRAN DATA LOAD SHEET III-5.1 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 Ti vs FT 60 
islooolhs6 Pal KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
aa cast PROG 


SAy) 


Pada Pa ETP CO OOO KORY OOOO 


Oe ea Meh Okt he et Be) 


o'e’s) Gove es 0 0 6 8 sO 8 + 6 6 0 0+ eee 
re MED SPITS A a eR OP re 
Ot ta ee Gap sve rerere---0.0 6-6 0 0.0 0-0 8 80 O08 
70.20.80. 8 + 0,0. 0,8," ee ee 8 oe 8 6 GD O.9 8 0 6 0.6.0 5 0.0.8 6 = 60.8 
RO Pali ee eM 
oe 


pe oe [ob Se 
Oo wy pu wy 
@ ee o~ 


ae 
oO 
es) 


ee 
oy 1m 


Pi 4617 
|. 4,6,8 
1 14,6,9! 
eR! 
eer 
oy 1447.2! 
Py 14,7,3) 
Ly 14,74 
wee 


MS OF ay IN THE FIRST TERM OF THE SPECTRUM GENERAL EQUATION. 
VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 7. 


Ww 


>> 
| 
>) 
uw Ww 


39 


B 
™ 


Sry, ~ 
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FORTRAN DATA LOAD SHEET III-5.2 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 TI 7S. 77 60 


000 KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


quan | toc jt] value [HE 
Lr 


2 
b 


11 15114 


51,4 


Say, = RMS OF Ay IN THE SECOND TERM OF THE SPECTRUM GENERAL EQUATION. ENTER 
Sry. * Ony, IF SECOND TERM IS NOT USED. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 7. 
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FORTRAN DATA LOAD SHEET III-5.3 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 TI Ta 77 80 


KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 


R R CASE PROG 
| 


Segm 2] 
4 
| 6 
| 28 
29) 
| 30) 
ma a 
ae Pee 
P14, 5.48: 
Yd, 54,9) 
| 1, 5,50) 
NG, 5,541 
Py 5452 
ey: PC) 
a 39 
Ted 45,555 0 


Say, = RMS OF Ay IN THE THIRD TERM OF THE SPECTRUM GENERAL EQUATION. ENTER 
Says = ayy IF THIRD TERM IS NOT USED. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 7. 
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FORTRAN OATA LOAD SHEET III-9 Prepared by 
Date 


STANDARD DATA INPUT 1 


68 Oo 7! 73 «77 


KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 
Rem cast PROG 


a ee a 


KR = GUST RESPONSE FACTOR *dny/oy. A CONSTANT IN THE GUST EQUATION SPECTRUM 
INPUT. 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 9. 
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FORTRAN DATA LOAD SHEET V-1.1 Prepared by 
STANDARD DATA INPUT 1 Date 

69 70 TI 73. 77 60 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
R R 


S-N DATA TABLE 1; N= f(Y,X) = CYCLES TO FAILURE 
(Y,X)= LOADS, SEE DATA SHEET II-5 FOR LOADS FORMATS. 


Nuay = MAXIMUM N CONSIDERED FOR DAMAGE CALCULATION, CYCLES 
= NUMBER OF Y ENTRIES, 2 < ny < 15. 
NUMBER OF X ENTRIES, 4 < ny < 15. 


Y = S-N CURVE CONSTANT LOAD, ENTER Y VALUES IN ASCENDING ORDER. 
X = S-N CURVE VARIABLE LOAD, ENTER X VALUES IN ASCENDING ORDER. 


THE ABOVE VALUES AND THE CORRESPONDING N VALUES, TO BE ENTERED ON SUBSEQUENT 
LOAD SHEETS, MUST BE INPUT ONLY WHEN IA AND/OR 14 = 1, 7, 13 OR 19 IS USED. 


= 
n 
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FORTRAN DATA LOAD SHEET I-2 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 Ti TS V7 80 
| 500 1 KEYPUNCH: DO NOT PUNCH BLANK DATA 
RR CASE PROG FIELDS. 


REFERENCE RUN NUMBER. 
CASE NUMBER. 
IRR AND ICASE NEED TO BE ENTERED ONLY IF ANY OF THE OPTIONAL PRINTOUT 
FLAGS IW1, IW2, IW4, IWS = 1. 
IEND = NUMBER OF SEGMENTS IN THE CASE (1 < IEND < 40). 
14 = FLAG WHICH SPECIFIES WHETHER GAG DAMAGE IS TO BE CALCULATED. 
= 0, NO. 
= 1 TO 24, YES. THE NUMBER CORRESPONDS TO THE S-N DATA TABLE NUMBER TO 
BE USED IN GAG DAMAGE CALCULATION. 
KEND = NUMBER OF THE LAST SEGMENT TO BE USED IN THE GAG CYCLE SPECTRUM 
DEFINITION. ENTER ONLY IF I4 # O. 
NEND = NUMBER OF GAG CYCLES TO BE DEFINED FROM THE INPUT SPECTRUM. ENTER ONLY 
IF 14 # 0. 
1 as STRUCTURAL ELEMENT ULTIMATE STATIC STRENGTH. 


S = WING AREA, Fi“. ENTER ONLY IF M3 = 8 OR 13 IN ANY SEGMENT. 
t = WING MEAN AERODYNAMIC CHORD, FT. ENTER ONLY IF M3 = 13 IN ANY SEGMENT. 


IRR 
ICASE 


S 


129 


FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page____—ioff 


FORTRAN DATA LOAD SHEET II-1 Prepared by 
STANDARD DATA INPUT 1 Date 
69 7O TI 72 TF 60 
[1 5]001 “ha KEYPUNCH: D0 NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 


Dy \/ 
LOG |» @ 1a pile ae 
1.6.97 Dix 
XX 
3] XA 
ma 
ey 


M3 = FLAG WHICH SELECTS THE SPECTRUM CYCLE INPUT FORMAT; A VALUE (1 TO 15) 
MUST BE ENTERED FOR EVERY SEGMENT BEING ANALYZED. 


= 1706: cn TABLES 1 TO 6. 
=o 2 
= 7: GENERAL EQUATION, rn = (IN or (ay)?/2lony)¥y7, 4 we 1, 2, 3 


= 8; GUST EQUATION, In = (IN ~Ay/b4Ay + i= 1, 2; A IS CALCULATED 
BY THE PROGRAM FOR et anhe. G. ver TioAL Gust LOAD FACTOR, 


= 9; GUST EQUATION AS ABOVE; A IS DIRECTLY INPUT. 

= 10 TO 12; (Sway SMIN? n) TABLES 1 TO 3. 

= 13; GUST EQUATION AS ABOVE; K,,, AND A ARE CALCULATED BY THE PROGRAM FOR 
AIRPLANE C,G. VERTICAL GUST LOAD FACTOR, 


= 14; GUST EQUATION AS ABOVE; A IS CALCULATED BY THE PROGRAM 
FOR AIRPLANE C.G, SIDE GUST LOAD FACTOR, 


= 15; GUST EQUATION AS ABOVE; A IS CALCULATED BY THE PROGRAM 
FOR VERTICAL TAIL SIDE LOAD. 
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FORTRAN DATA LOAD SHEET III-1.2 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 TI 73 77 80 
KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 
_F Pee i a 
Segm 3,478 4 ce re 
| 3,4,7.9 it oer ae © 
|g ,3,4,80 It} { 23] 3500 169.270 Ty 
Lol BET Id} [ 2af 3,501] 18,8,7,7 Fz 
| 3,4,82/ 16415 iy} {asf 3so2t Pay 
ee Ce ee eee ie ul 
ae CC LS es rk : 
|g: 3,4,8.6) 48 Ss eee 
ae Ee Ce ee ee i 
eet PLE ¥ 
a EI) 4,89 [5 
| 13 .3,4,9.0] 16828 +p 
1g 3.4,91] 16828 [| 
ma” ee a 
| 14 3,4,9,3] 189.27 aa 
ie Pe oe 
ae ees Pees et E 
me Re eer ii 
| ag 3.4.9.7] 189 27 > | eR 


F = A MULTIPLICATION FACTOR IN THE CYCLIC LOAD EQUATIONS. 
A VALUE MUST BE ENTERED FOR EVERY SEGMENT BEING UTILIZED, EXCEPT FOR 
SEGMENTS WITH M3 = 10, 11, OR 12, ONLY WHEN L1 = 1. 
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FORTRAN DATA LOAD SHEET III-1.3 Prepared by 
STANDARD DATA INPUT 1 Date 

69 70 TI 73 «77 60 

KEYPUNCH: D0 NOT PUNCH BLANK DATA FIELDS. 


Quan | Loc |#] VALUE 


? 

| 3 61.9 .6,2,5 
| 3,6,2,0/41,6,2,5 
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3 | 1316,2,5 
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P = A CONSTANT IN THE CYCLE LOADS EQUATIONS. 
A VALUE MUST BE ENTERED FOR EVERY SEGMENT BEING UTILIZED, EXCEPT FOR 
SEGMENTS WITH M3 = 10, 11, OR 12, ONLY WHEN IW3 = 1. 
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FORTRAN DATA LOAD SHEET III-1.4 Prepared by 
STANDARD DATA INPUT 1 Date ae 
69 70 Ti ge a 60 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


taetete 
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28 
| 29 


— 
My 


ee go 
rh 
—— aS 


=c 
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><, 
BS 


20 

° 

on 
Sots 
>< >< 
B<B< 


i 

00 [oo 
o fo 
<0 [oo 
i 
Se 
Eco 


9 sya 
P20) 18,12 of 


N = NUMBER OF LOAD LEVELS ENTERED IN THE LOAD SPECTRUM FOR ONE SEGMENT. 
(2 < N < 25) WHEN M3 = 1709 0R 13 TOI5;(1<n <25) WHEN M3 = 10, 11 
OR 12. 
A VALUE MUST BE ENTERED FOR EVERY SEGMENT BEING ANALYZED. 


= 
><D< 
SBS 


WO 
— 


40) .4.6352 
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FORTRAN DATA LOAD SHEET III-1.5 
STANDARD DATA INPUT 1 


69 70 TI Ts 677 


EP TOTO ECR” FOIE a Be WAN 
PEPTIC MPT TI RR rh he eae CMa) ae Ke 
seeeveyretetetetatitetetetateracacerene’ Ge .srte Me sisrersrrststastatatateraterereseceesregereststsSatatgtscecetacarererec ss *ste'Berereterecereres 
wittratste ele MececersMertetetatatatatetatetatetateretaretscerersgeseenstsestyteetatetecenatertsDies's'.' Datctatecateners 
pe eititacetesereteteteasgreratgrs coe crane! Malctscshtatetetateretetete rece enegeerseneeestatatetetecatarsceceseneneseDb enter i reratetgtecetere 
“srindetngraratatetessretetereteze, sory: Qereten: PecerecersresezerscnatarecesccatetetetetetatererecesetegrstatecstararctetecBpecscordiraress ak cee", 
RRMA MRE Pete AIM Ma NOR MIG HN MERI Co SOC Y 


ve ee eee 


PORE 
Td, 5,78] 101334, , ., | 10,0) 
a ve | 10)0| 
[ed , 5,8 0| 12,2,3,.6,! 
Tg, 5,81] 12,7235, , ,, | 10 
ft, 5,82] 16,7043 
Tg, 5,84] 12,7,2,0,9,6, , , | 10,5 
Td, 585] O48Ns .  | 145 
a, 5,86] 13,048% . | 10 
1, 587] 190489, Lt 
1a, 5.88[ 196089, |, | | 
Li) oa 


SPECTRUM, tn, 
A VALUE MUST BE ENTERED FOR EVERY SEGMENT BEING UTILIZED, EXCEPT FOR 
SEGMENTS WITH M3 = 10, 11, OR 12. 


KEYPUNCH: 


Page of 
Prepared by 
Date 


DO NOT PUNCH BLANK DATA FIELDS. 


igus. 
to 6.58...) 10 

reo, 100,658... | OO, 
pal ost ane dy.) OO 
oats 6.08) 128% ay yf 160 
30] 6.05) Ue 8.6 
al 606 ANTS 2... , | 00 
32 eto 132, | 100° 
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sof 2 Pt 
—asPoeaoll. ss tt 
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eve | ee Oe 

ol, etsd!.. tall, 


MULTIPLICATION FACTOR. 
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FORTRAN DATA LOAD SHEET III-2.1 
STANDARD DATA INPUT 1 


69 70 Ti 73 


Las 60 


RR CASE 


ite 


Bit 
eve 
eit 
EOE 
a42 
ty 
10,0 
ar 
¥a 
Lt 
a 
Lf 
wa 
a 
ae 
uF 
ae 
ry 
ae 
lot 
| 
‘o- 
int 


INPUT SPECTRUM Ay — In TABLE NO. 1. 
A MAXIMUM OF 25 VALUES MAY BE ENTERED. 


KEYPUNCH : 


Page of 
Prepared by 
Date 


DO NOT PUNCH BLANK DATA FIELDS. 


4 
3 


_— 


AB 


S 
© 


; i 
= ' 


AT LEAST TWO VALUES MUST BE ENTERED. 


Ay VALUES MUST BE ENTERED ONLY IF MS = 1 IS USED. 
In VALUES MUST BE ENTERED ONLY IF M3 = 1 IS USED. 
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FORTRAN DATA LOAD SHEET III-2.2 
STANDARD DATA INPUT 1 


69 70 7! 73 #77 60 
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— Pe fe 
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Fas, 


INPUT SPECTRUM ay — rn TABLE NO. 2. 
A MAXIMUM OF 25 VALUES MAY BE ENTERED. 
Ay VALUES MUST BE ENTERED ONLY IF M5 


= 2 
Zn VALUES MUST BE ENTERED ONLY IF M3 = 2 


KEYPUNCH : 


Page of 
Prepared by 
Date 


DO NOT PUNCH BLANK DATA FIELDS. 


Ss 
WN 


— | | 
So 1S 1S 
oO ini Pp 


ee — 
ieee — 
= 1c 


— 
— 
mn 


1,1,93 
] 2,0 


LLL + x) 


AT LEAST TWO VALUES MUST BE ENTERED. 


IS: USED. 
IS USED. 
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FORTRAN DATA LOAD SHEET III-4.1 Prepared by 
STANDARD DATA INPUT 1 Date 

69 70 TI 73 «77 80 

KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


Quan | toc 


No, = fn AT ay = O FOR THE FIRST TERM OF THE GENERAL OR GUST EQUATION SPECTRUM 
In INPUT. 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 7 TO 9, OR 13TO15., 
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FORTRAN DATA LOAD SHEET III-4.2 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 TI 73 77 60 
‘i sloo! KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 


Se 
398) WS... | 183 
Leeann 


No, = <n AT Ay FOR THE SECOND TERM OF THE GENERAL OR GUST EQUATION SPECTRUM 
En INPUT. ENTER No, = 0 IF SECOND TERM IS NOT USED. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 7 T0 9, OR 13 TO1S5. 
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FORTRAN DATA LOAD SHEET ITI-6.2 Prepared by 


STANDARD DATA INPUT 1 Date 


69 70 7I TS V7 80 


KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 


RR CASE PROG 


+| VALU 


OOOO DODO OOOOOD INOS OOOO OOOOOOOOOUUUUE 
BO a ee 


Ce CO Ce) ee 
Oe ee ee es ae ie OO ec a a) 
. 9 

Oe ee eee ed arene a ae ee eae 


penne” 
oe" “s 
oe ce 
or etere oe 
oes 0’ e's 
ee’ e's 
os 

. 
See 
eee "e 
eee 
-"* . 
ee ee 
oes" e'e*s 


1 16,82} 16,5,7 


oO» 
ee) 
nn 
Oy 
vy 
a | 


i 
7 


m = WING LIFT CURVE SLOPE PER RADIAN. 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8, OR 13. 
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FORTRAN DATA LOAD SHEET III-6.3 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 71 te te 60 
1 slooi{|l 6 PA: KEYPUNCH: D0 NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 


+] VALUE 


La, 69,7 


2.815 
2,78! 
2,60,4 
259.4 
23,13 


7 ,0,0 
: 7,0,1 
] F Oie 
7,0 
7 ,0 
7,0 


eae 
ee & Fe 


CO iN) [ow 


N N 
S ns 
W o> 


PK 
S jo fo 


co) el OS  E el a a aL. w 


~) 
co) 
eo) 


l250 745 


mn 


IATA FR 
— bs ie i ie 
Hr |W IN |e JO 
-Elsls 
| | 
if 


u 7 jl 


Ve = AIRPLANE SPEED, (EQUIVALENT AIRSPEED), KNOTS. 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8 OR 13. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA 
FORTRAN DATA LOAD SHEET III-6.4 


STANDARD DATA INPUT 1 


W 


Oe 


oe eee 7 
ee Pea Ee 
es oe Pees 6 
rerrirre Perri 
Ly fe) NSSF... | WOE! 
To | Serer i Pt 
fe hens oe! 
a eT eee i 
mS Terre ie 
1d, 5a 156213, ,, | 106 
nee fe 
a arr rh PPP 
| ag, 755] 1156213, ,, | 106) 


69 70 TI 73. 77 80 


KEYPUNCH: 


Page of 
Prepared by 
Date 


DO NOT PUNCH BLANK DATA FIELDS. 


PROG 


RR CASE 


AIRPLANE GROSS WEIGHT, LBS. 


oe | 
24! 
25, 
a) 
ai 
| 28 
29 
ee S01 «PRE 
3h 
| Se) y eT 
| 38! 9 ee 
| 341, 7.69 
Let a 
36 
a) 


=) 
Oo 
aN 


71. Taz 
38 5 shee 7,3 
40] 47,35 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8 OR 13. 
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FORTRAN DATA LOAD SHEET III-6.5 Prepared by 
STANDARD DATA INPUT 1 Date 

69 70 71 73: 77 60 

1 slo,01]f.6 Pal KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 


Rm CASE PROG 


PA 7st 19485, | 109 | 


Pa} aad t.... tty 
rT Tiaslt. tt! (as azastt 
oe) er 35 BO. 
Taal 5a... too) (shaadi. 
Ser TY Oe Te 
(ast hag ttt) [aztaursl Mt, | OO. 
PP neal toes TL aha 
Tid 1,7.22| 1326 joo] [a0] 1.74.2) 493 190. 
TT aaz2dt ta th <7€E 1 
ae CEP ee eee Ceres Lee 
Tag 7.251 2) 0] [oaa| nas. th 
ad dt] Cea 
edt. | casa 
Tad uz.za 322,,,,.../ 10] | selaurssl's yal 
Lf plete ey pe | Ty ee et | 
OCU Oe eg a oC Derren ee 
ad nat tt} (oe) azsall.. ph 


SIG = AIR DENSITY RATIO = p/p9- 
A VALUE MUST BE ENTERED ONLY FOR THOSE SEGMENTS WITH M3 = 13. 
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FORTRAN DATA LOAD SHEET III-6.6 Prepared by. 
STANDARD DATA INPUT 1 Date aan 
69 70 TI TS. 77? 60 
KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 


+/ € 


Lor L, Sas 


| Bal abba 
| _ Be 
26 

] 


S 
L 


ree | Pere 
Lererrrrr 


26 


See ieee 
291 3.66.6 

| 301 3.66.7 
Te 


. 
Ww 
oO 
@ 
> 
UW 
eo) 
5 


3,6,7,6 


( 


L = SCALE OF TURBULENCE, FT. 
A VALUE MUST BE ENTERED ONLY FOR THOSE SEGMENTS WITH M3 = 13. 


l = DISTANCE FROM AIRPLANE C.G. TO LIFT CENTER OF 
VERTICAL SURFACE, INCHES. 

A VALUE MUST BE ENTERED ONLY FOR THOSE SEGMENTS 
WITH M3 = 14 or 15. 
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FORTRAN DATA LOAD SHEET III-7.1 Prepared by 
STANDARD DATA INPUT 1 Date 
63 70 71 eS. fF 80 
KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 
| quan | toc |+; vatue [tle +} value sé | 
Ob p eee) TP, Or 2 ie 
cae 6 rr ee ee ee ee ee : 4, @ fF tot Lhd 10,0 
I 10, 


bi det: Cillian oP ab atgcda. 


Da ei & & £ & |0,! 


P. = PORTION OF FLIGHT TIME (OR DISTANCE) SPENT IN NON-STORM TURBULENCE. 
THE VALUE OF P IN THE FIRST TERM OF GUST EQUATION SPECTRUM INPUT. 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8, 9 OR 13 TO 15, 
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FORTRAN DATA LOAD SHEET III-7.2 Prepared by _ 
STANDARD DATA INPUT 1 Date 
69 70 TI v2 Fe 80 
KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 
RA CASE PROG 


4.22| 1000095 


hassel 
114) 8,4! 


rs = PORTION OF FLIGHT TIME (OR DISTANCE) SPENT IN STORM TURBULENCE. THE 
VALUE OF P IN THE SECOND TERM OF GUST EQUATION SPECTRUM INPUT. ENTER 
P, = 0 IF SECOND TERM IS NOT USED. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8, 9 OR 13 TO 15, 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA 


FORTRAN DATA LOAD SHEET III-8.1 


Page___—s off 
Prepared by 


STANDARD DATA INPUT 1 Date 

69 70 TI v3; 77 80 

(151007 KEYPUNCH: D0 NOT PUNCH BLANK DATA FIELDS. 

RR Case PROG 
| QUAN | LOC '+] ~~ VALU +| QUAN | LOC jt] VALUE = it E | 
p by fe ee ee mae hr 
Segm W .1,5,1,3) | ||, | sSeqm 21) 15,33) 129.5 1 1O7r | 
me ere (laf) | 22] 153.4 los’ 3 
9 15,15) 10992, | 100 23 |0,0 
er IN ee |3,0,5 1D 
ae eee eee o¢ 
— es 26/ .1,53,8| loVy 2 100 
| 8 i 28} 11,54,0] 1054.3, 1... | 140} 
ae ee : 100 | 
10 2,21 | , | 30/ 1,5,4,2| 125 lof | 
— 24, | _ 190) 33/ .1543| losy 3 | 100 
NF yn 2g: 108,69, 1 | 10,01 iQ | 
ar ss Joes. | ed 334 115.451 | ~~ 
| 4 152i 52s. | 1B! 34! 1,546] | in | 

19 1,5.2,7' (O97 21, . | 100] 351 154.71 | i} 3 

16 1,5,2,3: 14,68 | 10,44 361 ,1,5,4,8] | Be i 
PA a5, 2,01 1207 1G} {371 .145,4,9] | | tf 
|S gt i a! 128 | | its 4 
Fs se ai 13,2 a Set ee a we ayy | i | 
2 5,321 13.25,,,,,,; 160 101 waewl |» « « ¢-vmacn. | 


oe 
~_ 


NON-STCeiM TURSULENCE, 
TION SPECTRUM INPUT, 


= SCALE PARAMETER IN PROBABILITY DISTRIBUTICN OF RMS GUST VELOCITY FOR 
THE VALUE CF b IN THE FIRST TERM OF GUST EQUA- 


A VALUE HUST BE CNTERED ONLY FOR THE SEGMENTS WITH M3 = 8, 9 OR 13 TO 15, 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page____—off 
FORTRAN DATA LCAD SHEET III-8.2 Prepared by 
STANDARD DATA INPUT 1 Date 


69 70 TI ve 6TP 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


CASE PROG 


ee 10,543. 16,0 
| 26] 57,8 123,01... | | 100] 
bc OE ALB Ee eo yay OD 


281 1115 840) 1Q,0 

7,94, ag 15,81) 123,11. | 100 
3187 4 

W646 

391 sasieal (SSids cual) 0 


331 .1,5,85] | 

| 34] is. 8a | err 
35] 15.871 | | | 
361 15,88 


: 37} 115,89] | eer. ee 
P18 15479) 1833 38} 115,90] | | g 
Wl oard) 628 39/ 159 1{ | | 

=] — 
mer Le eee io] e321 i 


b, = SCALE PARAMETER IN PROBABILITY DISTRIBUTION OF RMS GUST VELOCITY FOR 
STORM TURBULENCE. THE VALUE OF b IN THE SECOND TERM OF GUST EQUATION 
SPECTRUM INPUT. ENTER bz = by IF SECOND TERM IS NOT USED. 


A VALUE UST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8, 9 OR 13 TO LS, 
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FORTRAN DATA LOAD SHEET I-2 Prepared by 
STANDARD DATA INPUT 1 hate 
639 70 7) pg: ee a 80 
KEYPUNCH: DO NOT PUNCH BLANK DATA 
RR CASE PROG FIELDS. 


IRR = REFERENCE RUN NUMBER. 

ICASE = CASE NUMBER. 

IRR AND ICASE NEED TO BE ENTERED ONLY IF ANY OF THE OPTIONAL PRINTOUT 

FLAGS IW], IW2, IW4, IWS = 1. 

NUMBER OF SEGMENTS IN THE CASE (1 < IEND < 40). 

FLAG WHICH SPECIFIES WHETHER GAG DAMAGE IS TO BE CALCULATED. 

0, NO. 

1 TO 24, YES. THE NUMBER CORRESPONDS TO THE S-N DATA TABLE NUMBER TO 

BE USED IN GAG DAMAGE CALCULATION. 

KEND = NUMBER OF THE LAST SEGMENT TO BE USED IN THE GAG CYCLE SPECTRUM 
DEFINITION. ENTER ONLY IF 14 # O. 

NEND = NUMBER OF GAG CYCLES TO BE DEFINED FROM THE INPUT SPECTRUM. ENTER ONLY 
IF 14 # 0. 

> = STRUCTURAL ELEMENT ULTIMATE STATIC STRENGTH. 


S = WING AREA, FT©. ENTER ONLY IF M3 = 8 OR 13 IN ANY SEGMENT. 
2 = WING MEAN AERODYNAMIC CHORD, FT. ENTER ONLY IF M3 = 13 IN ANY SEGMENT. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA . = 


FORTRAN DATA LOAD SHEET II-1 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 TI 73 «77 60 
1 s[O02 KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 


‘ 

Oo To Ie 
OO }|.o 

co fy 


o 
c 


27] 16,99 
| 2a 1.7.0.0) 
| 29f 3.7.0, 
80) oe 
| 31} 27,03 


WwW 
y ~ 


WwW) ww 


Ww 


FLAG WHICH SELECTS THE SPECTRUM CYCLE INPUT FORMAT; A VALUE (1 TO 15) 
MUST BE ENTERED FOR EVERY SEGMENT BEING ANALYZED. 


= 1 T0 6; tn TABLES 1 0 6. ; 
= 75 GENERAL EQUATION, tn = (1No,e7(4Y)*/2(oay)¥)7, 4 = 1, 2, 3 


= 8; GUST EQUATION, rn = (rN, ,Pye74¥/DHA)y gw 1 2. A TS CALCULATED 
BY THE PROGRAM FOR AIRPLAN# &.G. VERTICALGUST LOAD FACTOR. 


= 9; GUST EQUATION AS ABOVE; A IS DIRECTLY INPUT. 

= 10 T0125 (Suave Swry» m) TABLES 1 TO 3. 

= 13; GUST EQUATION AS ABOVE; K,, AND A ARE CALCULATED BY THE PROGRAM FOR 
AIRPLANE ©C,.G. VERTICAL GUST LOAD FACTOR, 


= 14; GUST EQUATION AS ABOVE; A IS CALCULATED BY THE PROGRAM 
FOR AIRPLANE C.G, SIDE GUST LOAD FACTOR, 


= 15; GUST EQUATION AS ABOVE; A IS CALCULATED BY THE PROGRAM 
FOR VERTICAL TAIL SIDE LOAD, 
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FORTRAN DATA LOAD SHEET III-1.2 Prepared by 
STANDARD DATA INPUT 1 Date 

69 70 Ti TS. tT 80 

KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 


RAR CASE PROG 


+| VALU 


| <tr — : Tj 
[Segm 1 Segm 21| ,3,49,8 11,2,.2,84, ., | 105, 
Td 3 479 We | aa ee SOA el LP 
|g 34,80 Yaa! Jaen eee | 1 
Tg 34,81] 112086, we eA yee Ld 
|g (3,4,8.2) 1,2,2,27 SPC arr ce 
[dg ,3,4,8,31 | Be BBO hy ole eee tL 
PLL Uae ee ble pe 

rh. | (2a) oS 508 Up aly py yee Ul) 

in ae ie eee 

Pe! 1 St Oe el eae 

COL) Lay eee | ie eed TP 

P| SO er ee oy | es 

= 

— Fo err 

tt ae eo eer re 

Pia) VS gp epee a 

rae Te | SeeeeetT nT ire 

ae Po} een ie 

10,5 oe errr ie 


F = A MULTIPLICATION FACTOR IN THE CYCLIC LOAD EQUATIONS. 
A VALUE MUST BE ENTERED FOR EVERY SEGMENT BEING UTILIZED, EXCEPT FOR 
SEGMENTS WITH M3 = 10, 11, OR 12, ONLY WHEN LI = 1 
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FORTRAN DATA LOAD SHEET III-1.3 Prepared by 
STANDARD DATA INPUT 1 Date 
69 7071 73 77 80 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


RR CASE PROG 


= 


-I184o | 
3,6,2,0\-118 40, 
ae eee 

i | 3,6,.2,3/711 4 9 
[ee 
sd 


100 
eta 
a 
as 
as 
tt 
ix 
ae 
if 
re 


36,0 9-I! 87,3 
) 3610-12228 


P = A CONSTANT IN THE CYCLE LOADS EQUATIONS. 
A VALUE MUST BE ENTERED FOR EVERY SEGMENT BEING UTILIZED, EXCEPT FOR 
SEGMENTS WITH M3 = 10, 11, OR 12, ONLY WHEN IW3 = 1, 
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FORTRAN DATA LOAD SHEET III-1.4 Prepared by 
STANDARD DATA INPUT 1 Date _ 
69 70 Ti te v7 60 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


— 


— 
ES CO 


it 
te 


& 


Lares i+ 


N = NUMBER OF LOAD LEVELS ENTERED IN THE LOAD SPECTRUM FOR ONE SEGMENT. 
(2 < N < 25) WHEN M3 = 1 TO 9 OR 13 TO15;(1 <n <25) WHENM3 = 10, 11 
OR td. 
A VALUE MUST BE ENTERED FOR EVERY SEGMENT BEING ANALYZED. 
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FORTRAN DATA LOAD SHEET III-1.5 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 Ti v3) 77 80 
[i slooz}f6 pal KEYPUNCH: 0 NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 


Loc jt 


VALUE 


iv 


T Fee : 
rl 
riot] {22}. 5.97 
7 5.78| 102.43, .- 1100] [23] 3.98 
Td szgf ozs | oe] [eal 5.99 
ds seo) 2933... | lox) | 281, 600 
Td ssf rout... ties} [ 26f. 60 lowes. 
on, 5,82] 05795, , 1195} | 27], woz 17.2,2,1,6,,,,] 04 
Td ss3fee7e2u 1105] | 26|, 6.03) 1aSc.c, ++] 100” 
9 sal 5312... 1105] | 20, 604) 1as,6e. TTL. 
Td 5,85] 155,1,1,07,,,1 105] [ 30] 605llo,,,.,. 1 Wb 
ut, 5.86 oa. tloo] (af, soll... .,.. [ab 
Ta 5.87] 04s, 1160] (a2), eo ig 
1s, 5,88] 21 5sa35,. [lo] [a3], soalliare 
ad, 589] sua... |p] [ sl, solane 
1, 5,99] 1543.,,,,,1100] | 35]. eaoll,,., 
a PA Oe ee eee 
T1789] SHsee, 4, 1165] al ett 
ad 5,93) 912,563, | 10,5 Leer 
ad, 5.95 ke 882, ,,,[lo4} [ aol, isl... 


T = SPECTRUM, £n, MULTIPLICATION FACTOR. 
A VALUE MUST BE ENTERED FOR EVERY SEGMENT BEING UTILIZED, EXCEPT FOR 
SEGMENTS WITH M3 = 10, 11, OR 12. 
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FORTRAN DATA LOAD. SHEET III-2.1 Prepared by 
STANDARD DATA INPUT 1 Date 

69 70 Ti TS: TF 60 

KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


Ot ee 
eva’ e eee ee ee 
erat ete ete ete ee ee a + 6 oe se am oe 
Sees ee ee 


. 
SP tt he ee Oe) Se 
se eee Rt ee Ce Oe 
avs rt 


y= 


= 
~ 


‘ 
2113 a 2 
i 


er 
3 
4 
| os 
a 
ald 
= 
Pas 
isi 
6 
| ay 
[ie 
ig 
20 
a" 
23 
rT 
2s 


> 
i 
m 


INPUT SPECTRUM Ay — =n TABLE NO. 1. 

A MAXIMUM OF 25 VALUES MAY BE ENTERED. AT LEAST TWO VALUES MUST BE ENTERED. 
Ay VALUES MUST BE ENTERED ONLY IF M5 = 1 IS USED. 

En VALUES MUST BE ENTERED ONLY IF M3 = 1 IS USED. 
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60000 te 
erere eee eRe ee 88 As Pe OOrOOIOO 
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eee tote state's tere” 
Re, LY Sete eh hee is ee eka a) MD ta 


FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page of 


FORTRAN DATA LOAD SHEET III-2.2 Prepared by . 
STANDARD DATA INPUT 1 )) als 
69 70 Ti 73 «TT 60 
15 KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


+ 
> 
ol 
<< 
m 
N 


Se 
s Witty i— is 


rotate Tre eta tate ete ate’ ep ete ep ete ep ate eta etn ets atatatetstetetn’s 
7. . * or Se eK WO) OC OO Ok” Chet Shela SI, 


is : 
=~] < 
SIS IGE 
N |S & 


SLs 
o) Ke, 


Nm 
W |W IW 

O~ (U1 
— Te 
ch 
hm |= 


m— fa Te fe 
2 1 1a 1© 
Oo 1 [xB PO 


marr wee Penner oe Tre s 
a aad 
el aed ee tea 
ae ee poe ie 
ST aa Lr ee 
ee ee i. eres 
7 ad 
— ne 1 mer 
peg Fo Tt Fas hs ef 
bn Me ews san w |e eh errr | a 
ee ee) Coa oe 
CAM) Cae 


INPUT SPECTRUM ay — rn TABLE NO. 2. 


A MAXIMUM OF 25 VALUES MAY BE ENTERED. AT LEAST TWO VALUES MUST BE ENTERED. 
Ay VALUES MUST BE ENTERED ONLY IF M5 = 2 IS USED. 
Zn VALUES MUST BE ENTERED ONLY IF M3 = 2 JS USED. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA 
FORTRAN DATA LOAD SHEET III-6.2 
STANDARD DATA INPUT 1 


63 70 TI TS. tF 60 


RR CASE PROG 


STDOOOOOOOD WOOD’ OOGOCOO OOOO rie 
OO Le Sa te a a CMa che hema ae 
RA at Sh hea MR COI Cheetah i ic 
ee MRM State Me He he ih ee Ce 
LKARRKK KS MS ge MN NN eh RR Se 
I MD MO Jt 


sb 
od 
nO 


~~ 


[ 


RIG! : 
= wie : 
J | Ww v* 


o 
wn 


o 
oOo 
tN 


H, 16,714 
1 164745 


6 ,7,5 


m = WING LIFT CURVE SLOPE PER RADIAN. 


KEYPUNCH: 


Page of 
Prepared by 
Date 


00 NOT PUNCH BLANK DATA FIELDS. 


KAO OtatK SOOCOOOO III 
HN OOH HC a IR I a a Sin aa 
wee see 6 6 8688 etete cre cree ne 0 8 0 Ol S eens eee 
stone tere ots Ble 0 c's 610.0 08.6.9 8.0.0.0 0.0 0.0.9,0,0,',9 95 %e ten eee ae ee Pyrat tr 
aes Par atete terete ee 5150 G0. 0.8.0.¢ 0 0 0 0 0M 6.) 6 ms 2 0.0 8» «8.8 
See » ere es ee Ot +10 0 6 0.8.0.0 6 0 86.0.0 . 
oes BO oe e's ee ove sete 0 0 8 5 0 9 08 © ats 
Ned a a ED . 


oe 


Ww |& 
W j|Ph 


wW 
> 


oOo 


— 
| 36 
38 


3 6,9,4 


6,9,5 


5 
Ww jh. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8, OR 13. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page of 
FORTRAN DATA LOAD SHEET III-6.3 Prepared by 
STANDARD DATA INPUT 1 Date 


KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 


Were eee ane og a ne Gin es 0 Mh a'e a a 6 be 6 es te te ee Oe 8 8 ee ee Oe ee ee Mew ew es oe ee 
OO a ee 

OO Cee) 
CN es <r ee ee ae 


Ce Re Re ee ae Be 
eS Se a 

Se ee eC aCe a ae 
eM Cee ee) 


14o 


29,05 
2,8,7,5 


. 
on 
- 
on 
— — — . 

. 

.- 

. 
ee 
. 

. 

“te 
“- 
‘e*e 
- 

.* 
ore 
-"« 

o 
oe 
on 

. 
"e 

oe 

4 

. 

oe 

oo 

oe 

On 

es 

* 

** 
oe 
ae 

“* 

Oo 
5865 
ee 
oe 
<<. 
tee’ 
oe 
“ 
ae 
s 
=. 
e*« 
ee 
os 
oe 
6O 


j 
ee ee ee 


Ve = AIRPLANE SPEED, (EQUIVALENT AIRSPEED), KNOTS. 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8 OR 13. 
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FORTRAN DATA LOAD SHEET III~-6.4 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 TI To. Fe 80 
[islooz|fi6 PA KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 


= 

PE my 

a. 
a pert = 
= a > ed, a a ae ae Ee 


L% >) 
Ps | 
<= 
—~ 


3} 1182588 


2 
rrr 
Per Peer 
Pet 
rr 
RPT) 


-S 
ok 
6 
i mg 


11 79.7.64 


W = AIRPLANE GROSS WEIGHT, LBS. 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8 OR 13. 


( 


Ui 
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FORTRAN DATA LOAD SHEET II1I-6.5 Prepared by 
STANDARD DATA INPUT 1 Date 

69 70 TI 73 «77 80 

1,5]002. KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


RR Case PROG 


POO AOE OO OOOO OOO IO OOOO OOOGOOOGO 
ee ee he ee ee Cae ee ae 


ee eC ek et ee Oe Oe a a a Ce) ee 
a te a ete ee ee ee + sei e ee ee Me ees Be ees ote +s toe so) 8 6 6 6 80 6 8 8 0.8 8 6 oF 8 Fo 8 4 8 we Be + o +B +o 6 oe + 88, 
a ke ay 
a a —* 
ve On) Mo 
Cr Se « 
. om *. 


ae A el i= 
Pete’ 5 prereset ey ee ee eee 7 

bat | 24} 1i7,3,6 13,5) 0 | 10,0 | 

eae er Leererer i 

fie ee Leen ee 


ae eee ire 
me OP ere 
30] 742 78.4, | 09 
ae ere ie 
P32) 7s 
a gh Ey ee pees) hy | 
mr eee ee 
mT eee ie 
$6 1,7,4,8{ | 
37/1 ,1,7,4,9] | | 
38] ,1,7,5,0] | ic 
39] 1,7,5,1] | 
40] ,1,7,5,2| | ai 


SIG = AIR DENSITY RATIO = p/po. 
A VALUE MUST BE ENTERED ONLY FOR THOSE SEGMENTS WITH M3 = 13. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA 


FORTRAN DATA LOAD SHEET III-6.6 
STANDARD DATA INPUT 1 


69 70 TI 73. 77 80 


Case PROG 


'Seqm 1 
a: 


| 36,40 
L 4 36.41 
| A 3.64.2 
a: 
— , 


raj ~ 
O1Nn 


: 6, 4 


5 
ee er 
10 3,6,4,7| 125 


Ww 
Oo 
wae 
oO 


TX) 
)) 
G 
o 
oo 


ae 

cet 

rizsoo.. [1b 
ma 

co 

co 

n 

lO,Y4 


SCALE OF TURBULENCE, FT. 


KEYPUNCH: 


Page of 
Prepared by__ 
Date 


DO NOT PUNCH BLANK DATA FIELDS. 


ritaz0,,,,,] 
23 
24 
25 


NM 


. Ss 
PEEP BE 


o 
Ww 


hw, 


— 


| 23) 
24. 
| __2s' 
| 26 | 
call 
— 
| 3 
p32 


‘ 
wt 


A VALUE MUST BE ENTERED ONLY FOR THOSE SEGMENTS WITH M3 = 13. 


DISTANCE FROM AIRPLANE C,G, TO LIFT CENTER OF 
VERTICAL SURFACE, INCHES, 


A VALUE MUST BE ENTERED ONLY FOR THOSE SEGMENTS 


WITH M3 = 14 or 15. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page oF... 


FORTRAN DATA LOAD SHEET III-7.1 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 Ti 73 77 60 
KEYPUNCH: D0 NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 


VALUE VALUE séd'+ 


1 18,0,6 


one 


P= PORTION OF FLIGHT TIME (OR DISTANCE) SPENT IN NON-STORM TURBULENCE. 
THE VALUE OF P IN THE FIRST TERM OF GUST EQUATION SPECTRUM INPUT. 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8, 9 OR 13 TO 15, 


16] 


FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page____—sioff 
FORTRAN DATA LOAD SHEET III-7.2 Prepared by 
STANDARD DATA INPUT 1 Date | saan 


69 70 71 vs 7 60 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


149,5| ll |\O/ 


j sonore 149,6| 1002.9, | 10.0 
1,4.9,7 | lO! 
Ye 
A 1 4,8,0] 10,6,02.8 0101 | Or 
PF 144, 8,1] 10,0,0,1,1 Ce eee 
| -10,—51,4,8,2| 10,0,0075 1,5,0,2 
I 
1,4,8,5| |O0O 009 5 


| 

000,045 1,5,0,8 | 
7! ot S08 | 

38 | ,1,5);1,0 | 

39) 15.1.1 | 

| 


40} 15,14 | 


. = PORTION OF FLIGHT TIME (OR DISTANCE) SPENT IN STORM TURBULENCE. THE 
VALUE OF P IN THE SECOND TERM OF GUST EQUATION SPECTRUM INPUT. ENTER 
P, = 0 IF SECOND TERM IS NOT USED. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8, 9 OR 13 TO 15, 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Pages of_ 


FORTRAN DATA LOAD SHEET III-8.1 Prepared by 
STANDARD DATA INPUT 1 Date 
69 70 TI 73: Fe 80 
/s[002 KEYPUNCH: 00 NOT PUNCH BLANK DATA FIELDS. 
RR Case PROG 
| QUAN | LOC '+} VALU it1€ | | QUAN | LOC /+| — VALUE thE! 
>) fo es ee 
Segm 115.1 Seam 21] 1.5.33] |29 ry 
a Fe | mS ernie 
15.15) 4 | | 23} 1535] 1055.09 190} 


15,16! !06,157 ,,,11Q0 4! 153,61 13.05, . 11 | 


71243 io/| [sf asadt......./ 


i hs E 
| 26} 153,98) 1O5509, | | G0! 
1,5,1,0 !32 | 107 27! 153.9 lzes 101 | 
BL 11)5,4,01 | re 
29 | 11541) 10,5,5,0.4,,, 41 10) 


iiaemrin| min cuore 


LLL | LEE 


i arrears rs 

52:4 10,64 5.7, 1 1: 100 0,0 | 
13 gs 38, 1G. . I 1Ot 33 4,5| | ia. | 
he EeCrcir | ic fie 4 34! 1.5.4.6 | fi. | 
| 15} 1.5.27: 10498 i 100 | 35/ 15471 | it. | 
Ml gtr ee, Wy 90 ie | SON 36| 1,548] | : 4 
ae ree Ea iia Ie  eaerrern 
18 1,5, 3,01 12,59 L1Q1 rc errr Te 
1g 5.944 18 | 100 39/ 1,5,5,1| | | be a 
| ad .3.5,3.0} 13.25 lo 40} ss2tt.. sia fd! 


- 


, = SCALE PARAMETER IN PROBABILITY DISTRIBUTION OF RMS GUST VELOCITY FOR 
NON-STORH TURBULENCE. THE VALUE OF b IN THE FIRST TERM OF GUST EQUA- 
TION SPECTRUM INPUT. 
A VALUE HUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8, 9 OR 13 TO 15, 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA 
FORTRAN DATA LOAD SHEET II1-8.2 
STANDARD DATA INPUT 1 


Page of 
Prepared by 
Date 


469 7071 Ae a 80 
KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RAR CASE PROG 
an 7 va 
Ss Ie ye Hae ee ihe 
| 10, / 
| 
5 mm Z3 
STS 4 ed a5] 187,72 | Tere Ls ss 
157.81 i AOS Y | 0,0} 
10) oF) 1 670 Ist | 10 | 
Or} {28} 5.0 | 14 | 
|0.1 79 | 21,5 811 11,724 10,0 } 
| |O,| | 30) 1.5.82| IS , ria 
itor] {an L583] SETer TTT ele 
| 10,0 , ae 5)a4) 11,9,0,5.4, ., | | 10,0 | 
» 5 ij lot] { 33] 15.a5 ee sD 
Y eee errrrert ire te Leewewerrs © 
: 10,0 351 15.87) | Fe 
eT Cee CP 
V7 41151619 IA pp | IOI F371 4158.91 | rit re 
A 1,701 33; , ) | {0,/ 38} .1,5,9,0} | ss tt 
19 15:71 O22, 10,1 39/ 159, 111, | 1c 
' : —_— 
og «.1.5,7.2] 17.7.5 oe eed OG) 9S hn ow we 3 | 
b, = SCALE PARAMETER IN PROBABILITY DISTRIBUTION OF PMS GUST VELOCITY FOR 


THE VALUE OF b 
ENTER b> = by IF 


STORM TURBULENCE. 
SPECTRUM INPUT. 


A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS 
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IN THE SECOND TEPM OF GUST EGUATIGS 
SECOND TERM IS NCT USc0. 
1S. 


WITH M3 = 8, 9 OR 13 10! 
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CASE NO. 
INPUT 
S-ULT . 


be 


REFERENCE PUN NQ. 


CATA 


IS 


THE FOLLCWIAC DATA 


0. 


---NEN)--- 


C-BAR- 
17.87 


-—-— Ss--- 
1.099 2107.20 


0 
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TEND KEND 14 
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CASE NC. 


Ls 


REFERENCE RUN NO- 


THE 


1S INPUT DATA 
---N SUR 


DATA 
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CASE AO. 
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REFERENCE RUN NC. 
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CASE 
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REFERENCE RUN NOC, 
THE FOLLCWIN! 


IS IN®°UT ‘ATA 


DATA 
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REFERENCE RUN NC, 
THE FOLLCWING 
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APPENDIX A 
PROGRAM A6PA CHANGES 


This appendix outlines the changes made to program A6PA since its 
publication as program 16PA in Reference 1. The listing given in Section V 
reflects these changes. The changes consist of a minor internal program 
adjustment in segment spectra summation procedure and an addition to calculate 
gust airplane side load factor and vertical tail load spectra. 


Re Spectra Summation 


The need to calculate GAG cycle damge (14 # 0) in order to perform segment 
spectra summation (IW4 = 1) was eliminated. Segment spectra summation can be 
obtained even when GAG cycle damage is not calculated (14 = 0). 


2. Airplane Side Load Factor and Vertical Tail Load Spectra Calculation 


The calculation of airplane side load factor and vertical tail load 
spectra due to lateral gusts were added to the program. The calculations are 
performed and obtained in the following manner: 


Cumulative frequency spectrum is calculated as, 


8 -Ayj/b]A -Ay ;/boA 
En(Ay) j = No, Pie + Nop Poe T 


Airplane Side Load Factor, M3 = 14 


Ay = airplane side load factor 


[VeSwKa, / a9eu | | vey ' “Yer 


calculated by the program. 


= (88 (5-3 — 
Ir vee// ay 


(26/(o/e0)0 Ce ay $¢9) (K/14)° 
Pay 7] W = radius of gyration in yaw, in. 


> 
il i 


vas 
t 


as 
u 
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Ve = airplane speed, KEAS 


Sy = wing area, ft2 
W = airplane weight, Ibs. 
Cy, = side force coefficient, tail contribution; 
T 
per radian = 3, (C " Fy 
By aT 
o = fuselage bending factor due to tail lift 
Cy side force coefficient, tail off; per radian 
BTO 
o = air density at altitude, slugs/ft2 
Py = air density at sea level, slugs/ft@ 
C+ = mean geometric chord of vertical surface, ft. 
at = lift curve slope of vertical tail; per radian = CL. 
V 
= ¢ —_ 
YB do 
T {Sy (" “i | 

S+ = area of vertical tail, ft? 

= ° ° if 2 
law yaw moment of inertia, Ib-in 
1¢ = distance from airplane c.g. to lift center of 


vertical surface, in. 


Vertical Tail Side Load, M3 = 15 


-Ay;/b]A -Ay;/boA 
En(ayj) = (No, Pye + NO, Poe T 
where, 
Ay = vertical tail side load, Ibs. 
A = (v a+S 
= calculated by the program 
Ves ats St, Kg, = same as defined for M3 = 14 
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Input Data when M3 = 14 or 15 


The input data in general is the same as in the vertical gust Spectrum 
calculation M3 = 13. The main difference is in the data required to 
calculate A , which is: 


Load Sheet No. 


Ct 
St 


The affected input data load sheets (new data or notes) are presented 
on pages 24] thru 246. 
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When input data printout is requested, the new input data is 
printed under the following headers: 


Parameter Header 
Ct VT Chord 
St VT Area 
Ivaw IYAW 
“¥g, SFC-T 
c SFC-TO 
Y8r0 
at SLOPE 
Tt SCALE OF TURBULENCE 


The last two headers are the names of the other parameters (m and L) 
Sharing the same location numbers. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page____of 


FORTRAN DATA LOAD SHEET II-1 Prepared by 
STANDARD DATA INPUT 1 Date 
49 To 7) 73 «77 80 
KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RAR CASE PROG 


Loc |+! VALUE 


e fs 
29] 1,7,0,1 XX 
a Fey erent i+: 


cers Leen) +: 
7 ,0,4 


[3s 
[asl a zoel, | XX 
| ccassilll |XX 

eer f 
Rie wenawend 


FLAG WHICH SELECTS THE SPECTRUM CYCLE INPUT FORMAT; A VALUE (1 TO 15) 
MUST BE ENTERED FOR EVERY SEGMENT BEING ANALYZED. 


1 TO 6; mn TABLES 1 TO 6. ; 
7; GENERAL EQUATION, rn = (1Noe7(4Y)*/2Coay)¥)7, 4 = 1, 2, 3 


8; GUST EQUATION, rn = (xN,.P;eOY/PHA)T, «4 © 1, 23 A IS CALCULATED 
BY THE PROGRAM FOR AIRPLANE &.G. VERTICALGUST LOAD FACTOR. 


9; GUST EQUATION AS ABOVE; A IS DIRECTLY INPUT. 

10 TO 12; (Sway SMIN? n) TABLES 1 TO 3. 

13; GUST EQUATION AS ABOVE; Kou AND A ARE CALCULATED BY THE PROGRAM FOR 
AIRPLANE C,G. VERTICAL GUST LOAD FACTOR, 


14; GUST EQUATION AS ABOVE; A IS CALCULATED BY THE PROGRAM 
FOR AIRPLANE C,G. SIDE GUST LOAD FACTOR, 


15; GUST EQUATION AS ABOVE; A IS CALCULATED BY THE PROGRAM 
FOR VERTICAL TAIL SIDE LOAD, 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page of 


FORTRAN DATA LOAD SHEET III-6.2 Prepared bys 
STANDARD DATA INPUT 1 ee een 
69 Yo fi... FS 79 80 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


ao |O 
oy IO 
Cale os 
sls 
— —«(| E _—|O CO 
oOo 
ee) 
ees 


m = WING LIFT CURVE SLOPE PER RADIAN. 
A VALUE MUST BE ENTERED ONLY FOR THE SEGMENTS WITH M3 = 8, OR 13. 


a = LIFT CURVE SLOPE OF VERTICAL TAIL, PER RADIAN. 


A VALUE MUST BE ENTERED ONLY FOR SEGMENTS WITH 
M3 = 14 OR 15. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page of 


eee 


FORTRAN DATA LOAD SHEET III-6.6 Prepared by sy 
STANDARD DATA INPUT 1 Date = sit™” 
69 70 TI T3 FT 80 
KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
RR CASE PROG 


}  6,3,8 


I. = SCALE OF TURBULENCE, FT. 
A VALUE MUST BE ENTERED ONLY FOR THOSE SEGMENTS WITH M3 = 13. 


| = DISTANCE FROM AIRPLANE C,.G. TO LIFT CENTER OF 
VERTICAL SURFACE, INCHES. 
A VALUE MUST BE ENTERED ONLY FOR THOSE SEGMENTS 
WITH M3 = 14 or 15, 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page  ——— 


FORTRAN DATA LOAD SHEET III - 6.7 Prepared by 
STANDARD DATA INPUT 1 Date 
62 O FI 73 =T @0 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


; MOMENT OF INERTIA ABOUT YAW AXIS, LB-IN® 
hd A VALUE MUS1 BE ENTERED ONLY FOR THOSE 
SEGMENTS WITH M3 = 14 OR 15. 


MEAN GEOMETRIC CHORD OF VERTICAL SURFACE, FEET, 
ENTER ONLY IF M3 = 14 OR 15 IN ANY SEGMENT, 

AREA OF VERTICAL TAIL, FT. ENTER ONLY IF 

M3 = 14 OR 15 IN ANY SEGMENT. 


— 
iI 


CQ) 
cr 
il 


YN 
i 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page____—ioff , 


FORTRAN DATA LOAD SHEET III -6.8 Prepared by sy 
STANDARD DATA INPUT 1 Date 
69 70 Ti 73 677 80 


KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 


14418 


C = SIDE FORCE COEFFICIENT, TAIL CONTRIBUTION, PER RADIAN, 


A VALUE MUST BE ENTERED ONLY FOR SEGMENTS WITH M3 = 14. 
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FATIGUE DAMAGE CALCULATION PROGRAM 16PA Page____—iof 


FORTRAN DATA LOAD SHEET III - 6.9 Prepared by. _ 
STANDARD DATA INPUT 1 Date 
69 70 Tl Ty FT 80 
16 PA KEYPUNCH: DO NOT PUNCH BLANK DATA FIELDS. 
R R CASE PROG 


< 


g, 
ct 
O 
Ww 


| 


| 

| 

rr 

int 

7 

He 

et 29| 3,9,4,7 
cL “lsat... tix | 
re 

ia! 

rr 

ore 

1% 

ha 

ig | 


W 
Ne) 
as) 
ies) 


3.9.49 


mami oad 

sai IO WwW I 

Ww 

~O 

W 

— 
Co 
(ea) 
W 
~O 
mn 
> 


: 
FE 


39.3 6 
i i 
bn A ee ye gee ee | 39} 3,957 | 1 | 
3,9 3,8 CC Deere te 
C = Airplane side force coefficient, tail off, per radian, 
YBto A value must be entered only for segments with M3 = 14. 
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APPENDIX B 


PROGRAM A6PD ALTERNATE OPERATION 


To generate a spectrum loading sequence through program A6PD requires a 
unique set of A6PA and A6PD input data, see Section V. Section IX describes a 
method of entering this data for each individual spectrum through a keypunched 
card deck. However, if one desires to generate a large number of spectrum 
variations from a basic input data set, an alternate method of input data 
management, which is more convenient and flexible, is available through the 
use of the CDC utility, UPDATE. This appendix describes this approach. 


The CDC utility UPDATE allows the storing of the input data on one 
permanent file (disk). The data is to be divided into separate groupings, 
henceforth referred to as "decks". Through the use of UPDATE , any number 
of these decks can be retrieved to create a proper set of input data to 
generate one spectrum. 


Each group of data is given a deck name. A *DECK name card must be 
placed at the beginning of each data group before the entire set of data is 
input to UPDATE to create a permanent data file (NEWPL). When retrieving 
these decks, it is necessary to place the deck names on compile directive 
cards (*COMPILE) which are input to UPDATE . Use as many compile cards as 
necessary and separate the deck names with commas: 


*COMPILE name 1, name 2, name 5, etc. 


For program A6PD, in the first *COMPILE card only, precede the first deck 
name entry with TWOCC3. This states that no word description is used in 
A6PA data (data sheet I-1 in Reference 1) and no such data should be entered 
into UPDATE . The data will be read into program A6PD in the sequence of 
deck names listed on *COMPILE cards. Consequently, first list the A6PA 
data and then the A6PD data. For further information on the CDC utility 
UPDATE see Reference 4. 


One method of dividing the program A6PD input data for storing on the 
UPDATE file, used in Reference 2 to generate a large number of spectrum 
variations, is: 
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A6PA data: Reference Run (Case 000) data 
Case data 
A6PD data. 


Do not include the last A6PA data card (1 in ccl and 99999 in cc69-73) with 
the A6PA data but include as the first card with the A6PD data. 


To illustrate the use of UPDATE to generate a spectrum loading sequence, 
the test case of Section X from Reference 2, will be used. In Reference 2, 
over one hundred spectra variations were generated from a basic set of data 
and additional data to produce the variations. The following deck names 
were given to the individual groups of data: 


A6PA data: 
Basic (original): RRXX , Ref. Run = XX, Case 000 


CXXYYY , Ref. Run = XX, Case = YYY 


Variation: RRXX/SZZ, Ref. Run = XX, Case 000, Spectrum = 22 
CXXY/SZZ, Ref. Run = XX, Case = Y, Spectrum = #2 


A6PD data: S22 » Spectrum = 22 


Each spectrum was given an identification number ZZ . However, if applicable, 
data deck of any spectrum can be used to generate any other spectrum. It must 
be remembered that the A6PA data is read into the program by a special sub- 
routine INPUT] which has the feature of replacing previously read data with 

a later entry. The deck setup (list of control cards) required to generate 

the test case (spectrum 1 or BS] in Reference 2) using Douglas Aircraft Company 
CDC computer, running under KRONOS, is given on the following page. 
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Job Control Cards to Generate Test Case Spectrum 
Spectrum #1, BS1) Using CDC UPDATE Utilit 


SER ,ALO1. 

BS1,MFL155000,T ,10  ,L_ | 

ACCOUNT, , . 

CID, Name, etc. 

REQUEST ,TAPE3 ,NT ,LB=KU,D=1600,VSN=RESERVE,T =—=—s—| 

ATTACH ,LGO=A6PDLGO,etc. 

ATTACH ,OLDPL=A6DATA6,etc. Note 1. 
UPDATE ,Q,D,8,C=DATA. 

COPYSBF ,DATA ,OUTPUT. Note 2. 
REWIND DATA. 

LGO ,DATA. 

/EOR 

*COMPILE TWOCC3,RR15,C15001 ,C151/S1 ,C15002 

*COMPILE C152/S1,S1 


/ EOF 


Note 1. A6DATA6 is the complete UPDATE data filename. 


Note 2. Use this card only if you wish to get a list of the input data 
(card images) for the spectrum being run. The list will 
appear in the printout before the program output. 
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